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1. Introduction 



Spatiotemporal pattern formation has been extensively investigated in a variety of chem- 
ical, biological, and physical systems J 1 ' Since the observation of the Turing pattern and the 
spiral wave pattern in chlorite-iodide-malonic-acid^ and Belousov-Zhabotinsky (BZ) ^ reac- 
tion, respectively, chemical reaction systems have attracted much attention on studying the 
pattern formation. In general, chemical systems can exhibit three types of properties: ex- 
citable, bistable, and Turing-Hopf. Some of these chemical reactions are light sensitive, such 
as the Ru{bpy)\ + catalyst BZ reaction. The light-sensitive feature of the media makes it pos- 
sible to control the spatiotemporal patterns by using certain strategy J 4 ' The control strategy 
can be classified to either external action such as the periodical forcing, I 5-11 ! or internal one 
such as the time-delayed-feedback J 12-16 ! The time-delayed-feedback has been used widely due 
to its adaptive feature. For example, in excitable system, the rigid rotation of spiral can be 
stabilized by changing the domain diameter of feedback control J 13 ' In bistable system, time- 
delayed-feedback can control of the nonequilibrium Ising-Bloch bifurcation, which realizes the 
transformation between the spiral and labyrinth patterns J 14 ' 

The systems of Turing-Hopf type have shown interesting spatiotemporal patterns, such as 
the hexagon, stripe, spiral and antispiral patterns. Some patterns are of desirable, such as 
the Turing hexagon grown in polystyrene film J 17 ' While some patterns are of avoid. Many 
methods have been attempted to control the Turing and Hopf patterns. For example, time- 
delayed-feedback with appropriate intensity has been used to suppress or induce the Turing 
patterns. ^ Recently, the systems near the Turing-Hopf co dimension- two phase space have 
shown many fascinating patterns, such as the oscillatory Turing patterns obtained in a Brus- 
selator model J 19-21 ! It is necessary to investigate the control of the spatiotemporal patterns 
near the Turing-Hopf codimension-two phase space. In this paper, we study the control of the 
transition between the Turing and antispiral patterns by using time-delayed-feedback strategy 
in a FitzHugh-Nagumo model. The effects of the time delay on the Turing and Hopf modes are 
analyzed by treating the time delay as perturbation. Numerical simulations show the flexibility 
of the time-delayed-feedback on controlling the transition between the Turing and antispiral 
patterns. We observe several dual-mode antispiral patterns and discuss their origins. 



2. Model 



This work is based on a delayed FitzHugh-Nagumo model: 

u t = u-u 3 -v + D u V 2 u + F, (1) 
v t = e(u - a x v - oq) + D v V 2 v + G, (2) 

where the time delay is applied with the forms: 

F = g u (u(t-r)-u(t)), (3) 
G = g v (v(t-r)-v(t)), (4) 

here, variables u and v represent the concentrations of the activator and inhibitor, and D u and 
D v denotes their diffusion coefficients, respectively The small value e characterizes the time 
scales of the two variables. The system described by Eqs.(l) and (2) can be either of Turing- 
Hopf, excitable, or bistable type. In this paper the parameter a\ is chosen such that the system 
is Turing-Hopf type. The parameter a determines the position of the uniform steady states 
on the nullclines of Eqs. (1) and (2). Unless otherwise specified, we choose a =0 such that the 
uniform steady states are uq—vq—Q. The parameters g u and g v are the feedback intensities of 
variables u and v, respectively, r is the delayed time. In order to study the effect of the time 
delay on the patterns of Turing-Hopf type, we first execute the linear stability analysis. In the 
analysis, we treat the delay as a perturbation by expanding the feedback terms Eqs. (3) and 
(4) as: 

u(t - 
v(t- 

So, we obtain: 

(1 + rg u )u t 
(1 + rg v )v t 

It can be seen that the terms of time is rescaled and determined by the delayed time and 
the feedback intensities. In the analysis, we consider the parameter ao=0. So uniform steady 
state of the system reads: uo—vq—0. Then, we perturb the uniform steady state with small 
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u - u 3 - v + D u V 2 u, 
e(u — a\v — a ) + D v V 2 v. 
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spatiotemporal perturbation (Su, 5v) ~ exp(\t + ikr), and obtain the following matrix equation 
for eigenvalues: 
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So we get the following quadratic equation for the eigenvalues: 

AX 2 - BX + C = 0, (9) 

where, 

A = (l + rg u )(l + Tg v ), (10) 
B = (1 + rg u )(e ai + k 2 D v ) - (1 + t 9v ){\ - k 2 D u ), (11) 
C = e-(l-k 2 D u )(ea 1 + k 2 D v ). (12) 
The dispersion relations X(k) of the system are defined by the roots of Eq. (9), 

Xi2 = b± vw^ac (13) 

The positive real parts of X(k) give rise to Hopf instability and Turing instability, which 
occur at k=0 and kj^O, respectively. In the absent of time delay, the stability diagram of the 
uniform steady state in the (ai, e) plane is shows in Fig.l. The Turing bifurcation line (solid 
line) and the Hopf bifurcation line (dash line) cross at one point. Applying time delay would 
affect the two bifurcation lines. Here, we take the parameters set as (a±, e) = (0.5, 2) as indicated 
by the open circle in Fig. 1, which is beyond the Turing bifurcation line and locates on the Hopf 
bifurcation line. Our emphasis is on controlling the pattern formation by using the strategy 
of time-delayed-feedback near the Turing-Hopf codimension-two phase space. In the following 
we will discuss the effect of time delay on the Hopf and Turing instabilities, respectively. Our 
analysis shows that the time delay affects both the Hopf and the Turing instability. 

The threshold condition of Hopf instability reduced from the linear stability analysis can 
be expressed as 

af = { } + T9v \ , (14) 
1 e(l + rg u y 1 ; 

and the corresponding frequency of oscillation is 
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Figure 1: Bifurcation diagram of the model system. Other parameters are: D u =0.01, D v =0.1, 
a =0.0, r=0.0, g u =g v =0.0. 

Eqs. (14) and (15) show that the strategy of time-delayed-feedback plays an important role 
on controlling the bifurcation line and the oscillatory frequency of Hopf instability. It can be 
seen from Eq.(14) that applying positive feedback to variable u (i.e., g u >0) and/or negative 
feedback to variable v (i.e., g v <0) is functionally equivalent to increasing e and/or decreasing 
ax, and vise versa. We want to mention that in a special case of g u =gv, the feedback does not 
affect the bifurcation line of Hopf instability any more. However, it still changes the oscillatory 
frequency of Hopf mode as is illustrated in Eq. (15). Fig. 2 (c) shows the dependence of 
the oscillatory frequency on the feedback intensities when applying the feedback with identical 
intensities g u =g v . The oscillatory frequency decreases dramatically with feedback intensities. 
The stronger the intensities of the negative feedback are, the faster the Hopf mode oscillates. 
In the case of g u ^gv, the dependence of the oscillatory frequency on the feedback intensity g u 
and g v are shown in Fig. 2 (a) and (b), respectively. The solid line part in Fig. 2 represents 
that the system is beyond the bifurcation line of Hopf instability. When the feedback g u is 
applied individually, the oscillatory frequency increases from zero and then decreases with g u 
as indicated by the solid line in Fig. 2 (a). When the feedback g v is applied individually, the 
oscillatory frequency decreases with g v as shown in Fig. 2 (b). Therefore, the time-delayed- 
feedback provides a way to change the oscillatory frequency of Hopf mode even under the 
condition that the intensity of Hopf mode remains unchanged. 

Fig. 3 shows the dependence of the wave vector of the most unstable Turing mode on the 





Figure 2: Dependence of the oscillatory frequency of Hopf mode on the feedback intensities 
in three cases: (a), g v =0; (b), g u =0; (c), g u =9v=9- Other parameters are: D u =0.01, D v —0.1, 
a 1= 0.5, a =0.0, £=2.0, r=0.1. 



feedback intensities. It can be seen that when the feedback g u (g v ) is individually applied to 
the system the wave vector increases (decreases) with the feedback intensity. This reveals that 
the wavelength of the Turing patterns can be controlled by using appropriate time-delayed- 
feedback. In a special case of g u =9vi the time delay do not affect the wave vector of the Turing 
mode any more as indicated by the dash line in Fig. 3. In this case, the time delay is equivalent 
to rescaling time as shown in Eqs. (7) and (8) and keeping the ratio of diffusion coefficients 
unchanged. Thus the wavelength of Turing pattern is unaffected. 
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Figure 3: Dependence of the wave vector of the most unstable Turing mode on the feedback 
intensities in three cases: solid line, g v =0; dash line, g u = 9v=9', dot line, g u =0. Other parameters 
are: D u =0.01, D v =0.1, ai=0.5, a =0.0, £=2.0, r=0.1. 




3. Two-dimensional numerical simulation 

In order to illustrate the above analytical results on controlling the pattern formation by 
using the strategy of time-delayed-feedback, we have carried out numerical simulation in two 
dimensions. The grid sizes are 200x200 which represents domain sizes 20x20 s.u., and the time 
step is At=0.01 t.u.. The boundary condition is zero flux boundary. In the above analysis, the 
time delay was treated as a perturbation, therefore the delayed time r should be a small value. 
The feedback intensities g u and g v can be of any values. Eqs. (7) and (8) have shown that the 
delayed time r and the feedback intensity g u (g v ) are incorporated together, which means that 
the effect of small delayed time and large feedback intensity, for example (r=0.1, g u =10.0), is 



equivalent to that of large delayed time and small feedback intensity (r=10.0, g u =0.1). In the 
following we only show the numerical results with small delayed time r=0.1. 




Figure 4: Transition from Turing to antispiral patterns with increasing g v in three cases: 
a =— 0.1, (a)-(d); a =0.0, (e)-(h); a =0.1, (i)-(l). The feedback intensities g v in (a)-(d), (e)- 
(h),(i)-(l) are 0.0, 5.0, 6.0, 14.0, respectively. Other parameters are: _D U =0.01, D v —0.1, ai=0.5, 
£ =2.0, r=0.1, g u =0.0. 

Fig. 4 shows three cases for controlling the transition between the Turing and antispiral 
patterns when the feedback g v is individually applied to the system. In the case of a =— 0.1, 
Fig. 4 (a)-(d) show the transition from the hexagon to antispiral patterns with increasing g v . 
When g v =0.0, it is a hexagon pattern as shown in Fig. 4 (a). This is because that the Turing 
mode is positive while the Hopf mode is negative as illustrated in Fig. 5. With increasing the 
feedback g v , the Hopf mode becomes stronger and stronger as shown in Fig. 5, which induces 
oscillating patterns such as the antispiral pattern. The Turing mode and the Hopf mode will 
compete, which results in the coexistence of the stationary Turing pattern and the oscillating 
patterns as shown in Fig. 4 (b). If the Hopf mode is dominant one the oscillating pattern 
swallows the stationary Turing pattern gradually, and if the Turing mode is prominent the 
stationary Turing pattern will occupy the whole domain finally. 

Fig. 4 (c) shows a dual-mode antispiral pattern which exhibits the competition between 
the Hopf and Turing modes near the Turing- Hopf co dimension-two phase space. The antispiral 
pattern originates from the negative dispersion du / d\k\<0 and the faster bulk oscillation oj >ook 
as indicated by the dash lines in Fig. 5. Here, we focus on the core region of the antispiral 
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Figure 5: Dispersion relation of the delayed model in four case: g v =0, 5, 6, 14. The solid (dash) 
lines represent the real (imaginary) part of the eigenvalue. The numbers indicated in the label 
represent the values of g v . Other parameters are: D u =0.01, D v =0.1, ai=0.5, ao=— 0.1, £=2.0, 
r=0.1, g u =0.0. 

pattern. As is well known, the tip of oscillating antispiral is a topological defect, and the 
amplitude near the tip is zero. However, in the present case of competition between the Hopf 
and Turing instabilities, the amplitude near the tip is not zero yet. The core of antispiral 
pattern contains stationary hexagonal spots originating from the Turing instability as shown 
in Fig. 4(c). The wavelength of the hexagon pattern near the core is about 1.1 s.u. which 
corresponds to the wave vector of the most unstable Turing mode k—5.7. The amplitude of 
the hexagon pattern is lower than that of the antispiral pattern as shown by the surface plot in 
Fig. 6. Far away from the core region, the hexagon pattern is damped by the Hopf oscillation. 
This is similar to the dual-mode spirals observed by Mau' 22 ' and Kepper' 23 ' , respectively. With 
increasing g v , the amplitude of the antispiral pattern increases gradually. The hexagon spots 
near the core is suppressed by the oscillation of large amplitude. The core of the antispiral 
pattern returns to a normal one as shown in Fig. 4 (d). 

In the case of a =0.0, Fig. 4 (e)-(h) show the transition from the stripe to antispiral 
patterns with increasing g v . This progress is similar to the aforesaid transition. The core of 
this dual-mode antispiral contains stripe pattern as shown in Fig. 4 (g). In the case of ao=0.1, 
Fig. 4 (i)-(l) show the transition from the honeycomb to antispiral patterns with increasing 
g v . The core of this dual-mode antispiral contains honeycomb pattern as shown in Fig. 4 (k). 
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Figure 6: Surface plot of the core of a dual-mode antispiral pattern. The amplitude of the 
hexagonal spots near the core is lower than that of the antispiral pattern. The parameters are 
the same as that in Fig. 4 (c). 

In a word, we can control the transition between the Turing and antispiral patterns by using 
time-delayed-feedback. 

The dual- mode antispiral is similar to the observation by Yuan et al in a CIMA reaction. t 24 l 
They called it Turing-Hopf mixed state and attributed its formation to the 3D effect of the 
reaction medium, where the Turing pattern and the antispiral pattern occur in different places 
in the third dimension. However, in our case the dual-mode antispiral pattern occurs during 
the transition as increasing g v . In this progress the Hopf mode becomes stronger and stronger, 
and the amplitude of the antispiral pattern increases gradually. When the amplitudes of the 
Turing pattern and the antispiral pattern are comparable as indicated in Fig. 5, the Turing 
pattern sustains near the core of the antispiral pattern, which forms the dual-mode antispiral 
pattern as shown in Fig. 4 (c), (g) and (k). Far away from the core region the Turing pattern 
is damped by the Hopf oscillation. With increasing g v continuously, the Hopf mode becomes 
strong enough. The amplitude of the antispiral pattern is much larger than that of the Turing 
pattern. The Turing patterns are damped nearly in the whole domain. Thus we observe the 
normal antispiral patterns as shown in Fig. 4 (d), (h) and (1). 

In the above simulations, we have shown the controllable transition between the Turing 
patterns and the antispiral pattern by increasing the feedback intensity g v . As illustrated in 
Eq. (14), the effect of increasing g v is equivalent to decreasing g u . The aforesaid transition 



from the Turing patterns (hexagon, stripe, and honeycomb) to the antispiral pattern can be 
reproduced by decreasing the feedback g u . This has been confirmed by extensive numerical 
simulations (not shown here). Therefore, the time-delayed- feedback provides a flexible strategy 
to control the Turing and antispiral patterns near the Turing-Hopf codimension-two phase 
space. 

4. Conclusions 

We have studied the transition between the Turing and antispiral patterns controlled by 
time-delayed-feedback in a delayed FitzHugh-Nagumo model. We treated the time delay as 
perturbation and analyzed the effect of the time delay on the Turing and Hopf instabilities 
near the Turing-Hopf codimension-two phase space. The time delay affects the wave vector 
of the Turing mode and the oscillation frequency of the Hopf mode. Analyzes show that the 
effect of small delayed time and large feedback intensity on the system is equivalent to that of 
large delayed time and small feedback intensity. The time-delayed-feedback provides a flexible 
way to control the pattern formation. The transition between the Turing patterns (hexagon, 
stripe, and honeycomb), the dual-mode antispiral, and the antispiral is numerically studied by 
applying appropriate feedback parameters. The dual-mode antispiral patterns originate from 
the competition between the Turing and Hopf modes. Our results could contribute to the 
control of pattern formations in the light-sensitive BZ or chlorine dioxide-iodine-malonic-acid 
chemical reactions. 
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